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In loop quantum cosmology the quantum dynamics is well understood. We approximate the full 
quantum dynamics in the infinite dimensional Ifilbert space by projecting it on a finite dimensional 
submanifold thereof, spanned by suitably chosen semiclassical states. This submanifold is isomor- 
phic with the classical phase space and the projected dynamical flow provides effective equations 
incorporating the leading quantum corrections to the classical equations of motion. Numerical work 
has been done using quantum states which are semiclassical at late times. These states follow the 
classical trajectory until the density is on the order of 1% of the Planck density then deviate strongly 
from the classical trajectory. The effective equations we obtain reproduce this behavior to surprising 
accuracy. 

PACS numbers: 04.60.Kz,04.60Pp,98.80Qc,03.65.Sq 



I. INTRODUCTION 

Loop Quantum Gravity (see [1-3] for reviews) is a non- 
perturbative approach to the problem of quantizing grav- 
ity. An open problem is that of the semiclassical limit, 
i.e. are there solutions to LQG which closely approxi- 
mate solutions to the classical Einstein's equations? Al- 
though this remains an open problem, concrete results 
can be achieved in the context of Loop Quantum Cos- 
mology (for recent review see [4,5]) in the spatially flat 
case. Loop Quantum Cosmology is a symmetry reduc- 
tion of LQG, i.e. a quantization of a symmetry reduced 
sector of general relativity. There, it is the case that one 
can indeed find semiclassical solutions in the quantum 
theory that closely approximate solutions to the classical 
Einstein equations, namely, the Friedmann equations at 
late times. Our goal is to to not just find semiclassical 
solutions but we would like to be able to go further and 
find effective equations incorporating the leading quan- 
tum corrections. 

It is a common misconception that canonically quan- 
tizing general relativity would just reproduce Einstein's 
equations without any modifications. The classical equa- 
tions are in fact modified by quantum corrections. Us- 
ing the geometric quantum mechanics framework [6| this 
has been previously shown for a dust-filled, spatially flat 
Friedmann universe [7| in the so-called /xg framework (see 
reference 3 in [^). In this paper we work in the so-called 
pL framework and show that this is also the case for a 
spatially flat Friedmann universe with a free scalar field. 

We can find effective equations arising from the ge- 
ometric quantum mechanics framework [6| ('effective 
equations' for short). These semiclassical states follow 
the classical trajectory until the scalar field density is on 
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the order of 1% of the Planck density where deviations 
from the classical trajectory start to occur. Then there 
are major deviations from the classical theory. However, 
states which are semiclassical at late times continue to re- 
main sharply peaked but on the trajectory given by the 
'effective equations'. Thus, we can approximate the full 
quantum dynamics in the (infinite dimensional) Hilbert 
space by a system of 'effective equations', incorporating 
the leading quantum corrections, on a finite dimensional 
submanifold thereof isomorphic to the classical phase 
space. Furthermore, we know from comparison with the 
numerical work that the leading corrections arc sufficient 
to capture essential features of the full quantum dynam- 
ics g. 

We begin by briefly summarizing the classical theory 
in ini We then provide an overview of the framework 
for effective theories that we use to obtain the 'effective 
equations' in IIIII Then in the quantum theory, IIVI we 
choose a family of candidate semiclassical states, verify 
that they are indeed sharply peaked, and compute the 
effective Hamiltonian constraint. We obtain the 'effective 
equations' in FVl and show that they are self-consistent to 
within our order of approximation. Finally, wc discuss 
how much this approximation can be trusted near the 
bounce point in the Appendix and conclude. 

Additionally, in this paper we use the choice of conven- 
tion: I'i = Gh, ciight = 1, K = SttG, and FRW parameter 
A: = since wc work in the spatially flat case. 



II. CLASSICAL THEORY 

In the classical theory, the phase space F is consists of 
pairs {A\,Ef) where A\ is an su{n) connection and the 
Ef its canonically conjugate field. Due to isotropy and 
homogeneity we have a fiducial background triad °e,° 
and connection "w^. Additionally, we may fix a fiducial 
cell with volume °V and restrict our calculations to that 



cell. There, we can write A^ and Ef as [^ 



E^ = P^ 



(1) 
(2) 



where °q is the determinant of the fiducial metric. Due to 
homogeneity and isotropy all the nontrivial information 
in A^ and Ef is contained in the variables c and p. The 
Hamiltonian constraint in these variables for a Friedmann 
universe with a free scalar field 0, in the regime p > is 
given by, 



C 



K7^ 



rO' + 



1P| 

2 p2 







where p^ is the conjugate momentum of the massless 
scalar field </). 

For convenience we switch to variables given by /3 = 
cj ^fp and y = p2 in which the Hamiltonian is given by 



K7^ 2 1/ 



(4) 



These new canonically conjugate variables are related to 
the old geometrodynamics variables via 



We can verify that pO|) and (|lip are equivalent to the 
Friedmann equations for a free scalar field when written 
in terms of ordinary ADM variables. Looking first at (fTT|) 
we see that it gives that, 

/3 = 7-. (14) 

a 

Putting in this expression in ([5]) we obtain, 
3 



K 2 a"* 



(15) 



C3'j which can be rewritten as the Friedmann's equations. 



3^7 = K 



2a^ a^' 



(16) 



or in terms of the Hubble parameter and the scalar field 
density, 



H'^ip 



(17) 



Now putting in the equation for P into ^ we can obtain 
the Raychaudhuri equation, 



(3 = 1- 
( 

V = a 



(5) 
(6) 



With the inclusion of a scalar field, the symplectic struc- 
ture now has an extra part. 



n 



K7 



dp A dV + d(l) /\ dp^, 



(7) 



where the scalar field is and it's conjugate momentum 
is p^. The Poisson bracket on the phase space is given 
by. 



{/,<?} = 



KJ f df dg 
dpdV 



dg df 
dpdV 



df dg dg df 

d4> p0 d(j) dp^ ' 



(8) 



where the phase space F consists of all possible points 
{/3, y, (J)tP4>}- The allowed pairs of points in F are those 
satisfying the Hamiltonian constraint 



C{P,V,<j>,P4,)=0. 



(9) 



Recall, that for a function that depends on the canon- 
ical variables, its time dependence is given by its Poisson 
bracket with the Hamiltonian. Thus the classical equa- 
tions of motion are given by: 



/3 = {/3,C} 
V = {p, C} 
<^ = {0,C} 
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v 



4 y2' 



(10) 

(11) 

(12) 
(13) 



3— = — 2k/9 
a 

Now ((T2|) and (fT3|) give us, respectively, 

<t>=-s:, 

P2 



and 



= const. 



(18) 



(19) 



(20) 



since C does not depend on for a free scalar field. These 
are just the equations of motion for a free scalar field. It is 
to the Hamiltonian constraint and these 4 equations, i.e. 
(Uni), dn]), dni), Uni), that we wish to find the corrections 
due to quantum gravity effects. 



III. FRAMEWORK FOR EFFECTIVE 
THEORIES 

Fortunately, there is the so-called geometric quantum 
mechanics framework which gives us a framework in 
which we can obtain these effective equations. In this 
section we briefly review the framework for effective the- 
ories. A more through review may be found in [7[ and 
[6|. This framework is especially suitable because it pro- 
vides a direct route to the effective equations from the 
Hamiltonian constraint without having to deparameter- 
ize the theory. That is, in background independent the- 
ories there is no canonical notion of time so one usually 
deparameterizes the theory by choosing one of the fields 
as an internal clock and then proceeds by considering 



the dynamics of the other fields with respect to the field 
serving as a clock. In particular, this model is depara- 
materized and analyzed in Q using the scalar field 4> as 
an internal clock. 

We take a brief detour and present the main idea be- 
hind the procedure for the simpler example of a particle 
moving in a potential along the real line. The proof that 
quantum mechanics has the correct semiclassical limit is 
usually based on an appeal to Ehrenfest's theorem, which 
is usually expressed in the form, for a particle in a po- 
tential, 



771- 



9*2 



/dV{x) 
\ dx 



(21) 



However, this holds for all states, not just semiclassi- 
cal ones. To recover classical equations for (x) what we 
would like to do is be able to pull the derivative outside 
of the expectation value. 



d^{x) dV{{x)) 



dt^ 



d{x) 



(22) 



Where in pulling the derivative outside, we have replaced 



— with QTpr, since after taking the expectation value 



dx 



the right hand side is now not a function of x but of 
(x) . Additionally, we would like to know if wc arc able 
to express the corrections to the right hand side of ([2^ 
in terms of a corrected potential, 



d^jx) 



d 



d(a 



[Vi{x))+SVi{x))] 



(23) 



That is, we are looking for an equation of motion satis- 
fied by the expectation values which includes quantum 
corrections to the classical equations of motion. A priori, 
it is not certain if we can do this. However, the geo- 
metrical formulation of quantum mechanics provides a 
framework in which we can answer this question. There- 
fore, we quickly review this framework. A more thorough 
review can be found in Q. 

In the quantum theory quantum states are represented 
by elements of a Hilbert space TC. More precisely, the 
quantum phase space consists of rays in Ti. The space Ti. 
can be made into a symplectic space via its inner product, 
which can be decomposed into its real and imaginary 
parts, 

{^\4') = ^G{^,i') + ^n{^,^), (24) 

where G gives a Riemannian metric and Q a symplec- 
tic form on Ti. Given a function / on the phase space, 
one can construct the Hamiltonian vector field associated 
with / by using the symplectic form. 



^if) 



f^^^Vb/ 



(25) 



It turns out that the flow generated on the phase space 
by the Hamiltonian vector field X?.^ corresponds to 



Schrodinger evolution in the quantum theory generated 
by/i. 

Taking expectation values of the operators correspond- 
ing to the canonical variables provides a natural projec- 
tion, TT, from states in the Hilbert space to points in the 
classical phase space. Thus, the Hilbert space, TC is nat- 
urally viewed as a fiber bundle over the classical phase 
space r. 



TT-.n-^T. 



(26) 



Because of the fiber bundle structure, the classical phase 
space can be viewed as a cross section of the Hilbert 
space. 

We can use this symplectic form to define a notion 
of horizontal vectors in Ti. Namely, vertical vectors are 
vectors whose components are in the directions in which 
the expectation values don't change and horizontal vec- 
tors are vectors orthogonal to the vertical vectors. It 
turns out that one has to use f2 rather than G to de- 
fine orthogonality. That is, two vectors, (f> and ip, are 
orthogonal if 



f](0,V) = O. 



(27) 



To summarize, fl gives us a notion of horizontal vec- 
tors which we use to construct horizontal sections of the 
Hilbert space. 

Because of the fiber bundle structure, any horizontal 
section can be identified with the classical phase space. 
At the kinematical level there is no natural section for 
us to choose. However, later when we consider dynamics 
in section |V] we can look for a natural section that is 
approximately preserved by the flow of the Hamiltonian 
constraint in a precise sense. 

A priori, we have no way of knowing whether such a 
section exists. However, if we can find such a section 
then the quantum dynamics on such a section can be 
expressed in terms of an effective Hamiltonian which is 
simply the expectation value of the quantum Hamilto- 
nian operator Q. The expectation value yields the clas- 
sical term as the leading term and gets corrections due 
to quantum effects in the subleading terms. This is the 
key idea behind our calculation. 

This can be done exactly for, e.g., the harmonic oscil- 
lator and approximately for several other physically in- 
teresting systems, including a dust-filled Friedmann uni- 
verse Q- The main result of this paper is proving that 
this can also be done approximately for a Friedmann uni- 
verse with a free scalar field and thus obtaining effective 
equations for this model. In the remainder of this work 
we will call these 'effective equations' because they are 
not effective equations in the traditional sense. These 
are on a different footing from traditional effective equa- 
tions since there is no precise definition of approximately 
horizontal sections. Nonetheless, we show that we can 
obtain these 'effective equations' in a well-controlled ap- 
proximation 



IV. QUANTUM THEORY 

Recall that in the full theory the elementary variables 
are the holonomies of the connection and the electric 
fluxes. Recall that in loop quantum cosmology there ex- 
ists no operator c [9| corresponding to the connection, 
and thus P also by extension. However, the holonomy 
operator, exp(^/2c) = exp{^^/Af3) does exist. Therefore 
in the quantum theory we work with the algebra gener- 
ated by exp(|\/A/?) and V 



A. Coherent State 

In the previous section, we saw that the taking of ex- 
pectation values provides a natural projection from the 
Hilbert space to the classical phase space. These expec- 
tation values can be taken in any state. However, for 
our investigation of the semiclassical limit it is natural 
to choose a semiclassical state. Simplest candidates are 
the Gaussian coherent states. They are natural because 
we can choose them to be sharply peaked at classical val- 
ues of the canonical variables, i.e. with small spread in 
both canonically conjugate variables. Let v and v' de- 
note the parameters proportional to the eigenvalues of 
the volume operator, related to the volume eigenvalues 



V and V via V 



3 ; 



(2p)f^^;andr = (^) 



3 It 



K " "^^^ " ~ y 6 ' K 
Then a Gaussian coherent state (V'| 



v' where 
with 



3V3V3 

Gaussian spreads e and e^ in the gravitational sector and 
scalar field sector of the state, peaked at some classical 
values /3', y , 0',P0 is given by, 

xe-54(P*-P*)'e-»"^'(p*-P*)(w;p^| 

o^^tpn{P4,){v;P4,\. (28) 



As usual, tp is defined on a lattice and the summation 
index v runs over the integers. It is important to note 
that the spread e is not constant but is a function of the 
phase space point. We will provide the reason and more 
details on the functional dependence in the next section. 
Recall that solutions to the constraints do not lie in the 
kinematical Hilbert space, but rather in its algebraic dual 
[lot ■ Physical states, such as the semiclassical state given 
in (pS)) should thus lie in the dual space and thus we write 
it as a so-called 'bra' state (V'l . At first it would seem that 
working in the dual space would be unmanageable but it 
does not pose computational difficulties. However, be- 
cause one fortunately has the "shadow state framework" 
to carry out calculations [ll|-p^. 

We will use {tp\ as our semiclassical state and calculate 
all of our expectation values in this state. The interpreta- 
tion of the basis ket is that the universe in the state \v;p^) 
has physical volume v in Planckian units and scalar field 



momentum p^. Our choice for (-01 as a Gaussian coherent 
state is not the most general but, rather, is the 'simplest' 
choice one can make to obtain these 'effective equations' 
in a late-time, large-volume approximation. For more 
general considerations it is possible to characterize the 
state in terms of all its 'moments' as in [l^. The higher 
moments acquire their own equations of motion which 
must be solved simultaneously with the equations of mo- 
tion for the canonical variables. 



B. Operators in the Quantum Theory 

In order to take expectation values, we need to define 
the operators in the quantum theory corresponding to the 
canonical variables. We now construct an approximate 
operator for /3 in terms of exponentiated /3 variables in 
order to verify that the state \1/ is indeed sharply peaked 
at classical value of /?'. Our point of departure is the 
classical expression 



1 



iVAp 



LVA/3 



(29) 



which is exact in the limit v A/3 -^ 0. However, our 
experience in the full theory tells us that in the quantum 
theory we should not be taking this limit to but to the 
area gap A — 4v37r7^p. So in our quantum theory we 
take 



A 



'A 



iyA/3 _ „-■^^/Ap 



(30) 



Thus, our operator /3a agrees approximately with the 
classical /3 in the regime \/A/3 <C 1. The choice for /3 is 
not unique but this is the simplest choice which is self 
adjoint (others have been considered in the literature, 
e.g. [iJl). For the rest of this paper we use this as our 
approximate /3 operator and drop the subscript A. Its 
action on our basis kets is given by, 



P\v-,P4>) 



1 



i\v + l;p^)~\v~l;p^)). (31) 



The action of the other operators is straightforward, V 
and p^ act by multiplication. 



V\v-p^) = (^^y '^v\v;p^)^ (32) 



\v;P(i>) = p^\v;i 



(33) 



and since we work in the p^ representation (/) acts by 
differentiation 



h d 

i 



(34) 



C. Restrictions on Parameters in Coherent State 

Before we move to computation we impose som.e physi- 
cally motivated restrictions on the parameters appearing 
in the coherent state (pS)) . The first pair of restrictions, 



%/A/3' < 1, 



(35) 
(36) 



corresponds to late times V ^ l^ and a ^ 1 respectively. 
Namely, that the scale factor be much larger than the 
Planck length and that the rate of change of the scale 
factor be much smaller than the speed of light. We will 
see later that ([55)) also holds well even at early times. 

The next pair of restrictions demands that the spreads 
in V and (3 be small, -^ <C 1 and -^ ^ 1, or equiva- 
lently: 



u'e> 1, 
e < VA/3'. 



(37) 
(38) 



The last pair of restrictions on parameters demands that 



the spreads in (f> and ^ 
or equivalently: 



are small. 



< 1 and 



P0 



«1, 

(39) 
(40) 



We use these physically motivated restrictions in our cal- 
culations in the remainder of this paper. We now return 
to showing that ip is sharply peaked. 



D. Verifying that tp is Sharply Peaked 

Now that we have a candidate semiclassical state, we 
can calculate the expectation values of the canonical vari- 
ables and verify that the state ip is sharply peaked at 
classical values of the canonical variables. 

Calculating {(3) we obtain. 



(/?> = 



sin(-VA/3') 



(41) 



Thus we must take e ^ 1 for (/3) to approximately agree 
with the classical /3. Similarly, 



0') = X 



1 



' cos(VA/3') 



Thus, 



A/?2 



2e2 

A ' 



(42) 



(43) 



where to get the last line we have used cos(v A/3') ~ 1 
and e ^ 1 . Now we obtain the expectation value and 
spread of V, 



{V) ^ v, 



(44) 



where to obtain the last line we have done Poisson re- 
summation on the sum over n. Similarly, 

AV' = ^K^y. (46) 



iV' 



v- 



(45) 



2e2 ^2 I 6 



Thus, up to the approximation used to arrive at (j43|) the 
product of uncertainties is the minimum possible. 



A(3AV 



T2' 



(47) 



We make a brief detour to show that we can satisfy 
both of the conditions (P7|) and ([55]) on the phase space 
if the width e of the coherent state is properly chosen 
as a function of the phase space point. Notice that ([37]) 
requires that e satisfy v'e ^ 1 and (|38p requires that 
e ^ a/A/?, but V^P ^ 1 so it appears that there is a 
tension between these 2 conditions. Nonetheless, if pj^ 
is large, which is reasonable for a universe that would 
increase to macroscopic size then we can choose e(w') 
such that both of the conditions ^- and -^ are satisfied. 
Looking at the relative uncertainties, we get 

3 

<1 



AV 
V 
A/3 
T 



1 



'p 



/8^Y 
^/2ey' K\ 6 J 

A^ 



<1 



(48) 



(49) 



The first condition shows that we need to choose e as a 



function of the phase phase variable V' s.t. e 



(for 



A ^ {^p-)'^K ■^) in order to satisfy the inequality. As 
for the 2nd condition we use the the classical form for 

2 2 

/32 = ^^-p, where p = -^ for a scalar field, to show that 
the inequality is equivalent to taking v^SttS <^ p^. For 
a universe that grows to macroscopic size such that p^ 
satisfies this then we can choose A in such a way to satisfy 
A < {^fK-\ e.g. A = 70, p0 = 5000. Therefore we 
can satisfy these two conditions simultaneously. 

Continuing, taking the expectation values of (f> and p^ 
yields, 



('A) 



and 



so 



2^^ 



A(/>: 



i^ 

V2- 



A0 



Thus, to have -^ ^ 1 we must have e^ ^ 



for 



and 



P4,, 



{P4>) ==P0> 



(pI) 



12 
P<t> 



1 

2d' 



(50) 

(51) 

(52) 
Similarly 

(53) 
(54) 



therefore we have for Ap^ 



The action of the these on our basis kets is given by, 



Ap0 



y2e. 



(55) 



Ap,/. 



Thus, to have — ^ ^ 1 we must take p^e^ ^ 1 . 

Notice that this coherent state saturates the Heisen- 
berg uncertainty bound as it should). 



A(t)Ap^ 



1 



(56) 



Therefore we have confirmed that the expectation val- 
ues of the operators corresponding to the canonical vari- 
ables agree with their classical values. Additionally, we 
have also shown that this coherent state saturates the 
Heisenberg bound, as it should. Therefore it is a natu- 
ral kinematical semiclassical state to work with in order 
to obtain these 'effective equations' from the geometric 
quantum mechanics framework. 



E. Expectation Value of the Hamiltonian 
Constraint Operator 

Having verified that ip is indeed sharply peaked around 
classical values of the canonical variables, we proceed to 
calculate the expectation value of the Hamiltonian con- 
straint operator. Recall that, if we can find an approx- 
imately horizontal section of the Hilbert space then the 
effective quantum dynamics on that section is generated 
by the effective Hamiltonian which is just the expecta- 
tion value of the Hamiltonian constraint operator. We 
calculate this expectation value now, leaving the proof 
that the section is approximately horizontal to a later 
section. 

Now notice that the self-adjoint constraint is given by 



^ " 167rG^^''""+2^'^ 



(57) 



where 



^ qra 



Ca, 



sin(/ic) 



24i sgn(/.t) 
87r73/i3;2 



-(f) --(f) 
,cos(f)fsin(^)]}sin(/2c) (58) 

-i®pl (59) 



^grav\'V\Pij)) 

C^\v\p^) 



f+{v)\v + A-p^) + h{v)\v;p^) 
+f^{v)\v^A-p^) (60) 

3 

^' ^ ^^ B{v)pI\v-p^) (61) 



2 \?>T^lll 



(62) 



where 



/+(«) = T-a 



f-iv) 
fo{v) 



16V 6 7 

/+(«-4) 



^|f-f 2|||w-M|-|t; + 3|| (63) 



(64) 
(65) 



and 



B{v)=\7;] K\v\\\v + 1\^ ^\v-l\^f (66) 



and 



P ^\v;P4>) 



Sirjll 



B{v)\v;p^) 



(67) 



Using these one can go ahead and compute the expec- 
tation value of the Hamiltonian constraint (C). Note 
that such a summation includes a summation over nega- 
tive V even though we are considering late times v ^ 1. 
However the contribution introduced from these terms is 
exponentially suppressed. For a proof of this see Q and 
the appendix in [9| . We can evaluate the sum by Poisson 
resummation and find an asymptotic expansion for the 
first term in the series obtaining the effective constraint 
to leading and subleading order. 



(C) = - 



WnGj^fl' 



2 ,-,'2' 



l + e~^'' (2sin2(\/A/3')-l 




v' 



(68) 



F. Equations of Motion 

Recall that to evaluate O for some operator O we sim- 
ply take the commutator between O and the Hamiltonian 



and divide by ih. We begin with (/3), 

(V'lV') 



(/3> 



(69) 



which can be computed using a similar expansion to that 
used to find the expectation value of the Haniiltonian 



constraint to yield, 



J 



0) 




4e * cos 



^yA/3'')+4e-^^osQ^/A/?' 



r 



ie ■3'=' cos ( -%/A/3' j 



(70) 



Now for {V) 

{1P\1P) 

Notice, V commutes with C^ so there is no contribution 
dependent on the scalar field and we just retain the con- 
tribution from the commutator with Cgrav which yields 



{V) 



ZV _4^2sin(2VA/?') 



7 2VA 

Now focus on (0) and (p,/,), 



(</') 



since 6 commutes with C, 



(V'lV') 



grav 



(^)-f^ + ^(^ 



Similarly, 



{'P<t>) 



(^lik 



ih V167rG^ 



b,Ca 



\[p^,c^])m 



(^'k'> 



(72) 

(73) 
(74) 

= 

(75) 



since p^ commutes with both Cgrav and C^. 

From these expressions we obtain the 'effective equa- 
tions'. 



V. EFFECTIVE EQUATIONS 

Recall that in the geometric quantization picture we 
take as our basic observables, the expectation values. 

^ = (/3> = -^^'~''^ sin(iVA/?') (76) 



V ^ {V) = v 

^ = (0) = 0' 

P<l> = (P*) = P'4> 



(77) 
(78) 
(79) 



Recall, that the coordinates on the classical phase space 
are the expectation values, thus these barred variables are 
the coordinates on the classical phase space. So we search 
for an effective description in terms of these variables. 

We thus invert these ([75)1 - ((7^ for the primed vari- 
ables and express the evolution equations in terms of 
the barred variables. We look at the first few terms of 
these asymptotic expansions to get the leading and next 
to leading behavior and apply the approximations listed 
above in llVCI Doing this we obtain the 'effective equa- 
tions of motion' (the main result of this paper). 



C 



K7^ \ 4 / K7^ A 



+|t[i + o(y-^y-V^)] 



(80) 



/?= |-^l-iA/32[-2/3^-HA/5^] 



1 



Pi 



^ y'l - ^A/32-^[l + 0{V-\V-\~^)]m 



* = f + 0(v-») 

Pd = 



(82) 

(83) 
(84) 



As discussed in the Appendix, Eqs ( (|80p and ([5^ imply 



that the Hubble parameter H 
effective Friedmann equation. 



H' 



where 



f" 



Pcrit 



p 



Pcrit 



0{e') 



satisfies an 



(85) 



K72A 



(86) 



which incorporates the leading quantum corrections. 
Somewhat surprisingly, as numerical simulations show, 
these corrections are already sufficient to correctly repro- 
duce the main features of full quantum dynamics. These 
equations have been used in the literature to make cer- 
tain phenomenological predictions (see, e.g., |15l|). 



One may ask, what is the error involved in making 
the approximations that led to equations ([5D|) - ([M|) . In 
these equations wc are keeping corrections of 0(A/3^) 
relative to the classical expressions but ignoring terms 
of 0{A'^f3^). How much of an error is one making by 
ignoring these terms? Let us compute, A/?^, 



A/32=A72^ = ^ 1-^ 



Pc 



Pc 



(87) 



Notice, this function starts off at at late times since 
p/ Pc ~ 0, it slowly increases reaches a maximum value 
of 1/4 at p — pc/2 and then again goes to at the point 
p/pc = 1- Therefore, at the worst we are making an error 
of about 6% in ignoring terms of the order A^/3^ and this 
occurs at p = Pc/'2. At other times, i.e. at late times and 
at times near the bounce point, this approximation is 
very good. 

Additionally, one can puUback the symplectic structure 
fl to our candidate horizontal section and analyze the 
dynamics generated by this effective Hamiltonian. Does 
it give the equations ([81]l - ((84|) ? One can puUback H. and 
verify that it is given by 



n 



1 



K7 



1 - iA/52 



zdl3 AdV + d(t)/\ dp^ 



(88) 



and that its associated Poisson bracket is 



{fi9} = -TT 



1^1 rr 

Of dg 

d(j}dp4, 


4 ^ \dpdV 
dg df 
dp4, d(j) 


dgdf 
dpdV 



(89) 



We can verify that our 'effective equations' are consis- 
tent by checking that the Poisson brackets in terms of 
the barred variables hold, to within our order of approx- 
imation. That is, (/3, V, 4>iP4>) are the coordinates on the 
horizontal section, F, and we wish to know whether they 
are preserved, i.e. is the vector (/?, y, ^,^0) tangent to 
r or off it? If it is not tangent to F then is it approxi- 
mately tangent? That is, are its components off F small 
relative to the tangential components? This is indeed the 
case, since to our order of approximation the equations 
of motion in terms of the barred variable hold. Indeed, it 
can be verified that the following equations hold in terms 
of the puUback of the symplectic structure to the barred 
variables. 



V 



{P.P.C] 



(90) 



This is a highly nontrivial consistency check on the for- 
malism as well as approximations used at intermediate 
steps. 



VI. CONCLUSION 

Thus, we have constructed kinematical quantum states 
that closely approximate solutions to the classical Ein- 
stein equations. Therefore, at least for the case of a 
Friedmann universe with a free scalar field in the con- 
text of loop quantum cosmology, it is the case that there 
exist suitable quantum states that closely approximate 
solutions to the classical Einstein equations and remain 
sharply peaked along a quantum corrected, semiclassical 
trajectory. 

Furthermore we note, that wc have shown that the 
classical Einstein's equations are not just reproduced 
identically in the quantum theory but they do indeed 
pick up corrections due to quantum effects. In this work, 
we have used the geometric quantum mechanics frame- 
work and approximated the full quantum dynamics in 
the (infinite dimensional) Hilbert space by a system of 
'effective equations', incorporating the leading quantum 
corrections, on a finite dimensional submanifold isomor- 
phic to the classical phase space. 
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VIII. APPENDIX: VALIDITY OF EFFECTIVE 
EQUATIONS NEAR THE BOUNCE 

It is instructive to write the 'effective equations' in 
standard form. Therefore, in this appendix we address 
the question: what is the form of the corrected Fried- 
mann and Raychaudhuri equations? 

We can obtain the corrected Friedmann equation by, 
working to next-to leading order, solving (|80)) for /3, sub- 
stituting into (|5^ . and then computing the Hubble pa- 
rameter H = - = rXr, 
a 3V ' 



H' 



up 

T 



1 



p 



Pcrit 



+ 0{e^) 



where 



Pcrit 



K72A 



(91) 



(92) 



Notice, in addition to recovering the classical ^p term 
we get a correction term that becomes important when 
the density p becomes very large. We see that the minus 
sign allows _ff , and hence a to be zero. Thus allowing the 



possibility of a bounce. Thus, even though these equa- 
tions hold at late times since we are working in the semi- 
classical regime, we already see that gravity is becoming 
repulsive and allowing the scale factor to bounce when 
we evolve backwards in time and approach the classical 
big bang singularity. 

On (j9ip we must make two comments. First, a priori 
one would not expect (PT|) to describe the correct dy- 
namics near the bounce point because the bounce point 
pI Pc = 1 lies outside of the regime of our approxima- 
tion. That is, the natural domain of applicability of 
these 'effective equations' is a late-time, large- volume ap- 
proximation and the point pj Pc = 1 lies well outside of 
that regime since the condition ([38|) fails badly at that 
point. However, numerical work (see reference 2 in [8|) 
has shown that the dynamics is indeed described well by 
(j9ip even at the bounce and hence this approximation 
and the results obtained in this work continue to be good 
even beyond their expected regime. Secondly, near the 
bounce point the term in parentheses in (j9ip is approach- 
ing 0, but since the approximation (j38p is breaking down 
there it is not clear that the O(e^) term is negligible there. 
Thus (|91|) appears to break down there, but the numeri- 
cal work has shown that it holds very well with negligible 
O(e^) corrections. As is not uncommon in physics, the 
effective theory is valid well beyond the domain for which 
it was constructed. 

So now we have the effective Friedmann equation in 



this model. It would be instructive, for completeness, to 
also obtain the conservation equation and the corrected 
Raychaudhuri equation for this effective theory. Recall- 
ing that p = p for a massless scalar field, classically these 
equations are 



P = -6-p, 
a 

3— = — 2k/9 
a 



(93) 
(94) 



Let us compute the corrected versions of these equa- 
tions. Using the modified Poisson bracket, we can com- 
pute p = {p,C}. Doing this we obtain precisely (P5)) . 
so the continuity equation is not modified. Similarly, we 
can take ([5^ and compute V = {V,^}, expressing it 
in terms of the scale factor a, and solving for 3-. One 
obtains 



3^ = -2.p(l-|^)+0(e^) 
a \ 2 Pc ' 



(95) 



Notice that classically this expression is always negative. 
However, now there is a correction of 0( — ). Notice, that 



the corrected Friedmann equation (|9ip tells us there is 
a bounce at p/pc = 1. At that density the corrected 
Raychaudhuri equation (|95p tells us that a is positive as 
it should be if there is to be a bounce. 
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